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Abstract 

The effective action for light mesons in the external uniform static 
electromagnetic fields was obtained on the basis of QCD string theory. 
We imply that in the presence of light quarks the area law of the 
Wilson loop integral is valid. The approximation of the Nambu-Goto 
straight-line string is used to simplify the problem. The Coulomb-like 
short-range contribution which goes from one-gluon exchange is also 
neglected. We do not take into account spin-orbital and spin-spin 
interactions of quarks and observe the p and K* mesons. The wave 
function of the meson ground state is the Airy function. Using the 
virial theorem we estimate the mean charge radii of mesons in terms 
of the string tension and the Airy function zero. On the basis of 
the perturbative theory, in the small external magnetic field we find 
the diamagnetic polarizabilities of p and K* mesons: I3p = —0.8 x 
10-"^ fm^ (3k* = -0.57 x 10"^ fm^ 



1 Introduction 

One of the important problems of particle physics is the confinement of 
quarks. There is progress in understanding of properties of mesons as a sys- 
tem in which quarks and antiquarks are connected by the relativistic string 
with a Nambu-Goto self- interaction [1]. This binding interaction becomes 
strong at large distances and therefore it is impossible to describe the phe- 
nomenon using the perturbative approach. There are some difficulties in 
evaluating meson characteristics in the general case of a complicated string 
configuration. Naturally, as a first step, we make some approximations and 
model assumptions to simplify the calculations. So here we consider the 
straight-line string as a simple configuration and quarks attached to the ends 
of the string. Such configurations were studied in [2] . In the present paper we 
investigate mesons in external, constant, and uniform electromagnetic fields 
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and use the path integral approach. It should be noted that in potential-like 
models [3-5] meson characteristics are described reasonably. But in these 
approaches there are some assumptions: (i) the relativistic invariance is only 
the approximate, and (ii) the constituent quark masses are used (but not 
current quark masses). 

The recent development of the QCD string approach [6-11] showed good 
results in describing heavy quarkonia, baryons and glueballs. The QCD string 
theory takes into account the main nonperturbative effects of strong interac- 
tions: chiral symmetry breaking (CSB) and the confinement of quarks. Chiral 
symmetry breaking gives a nonzero quark condensate {{qq))- As a result the 
light quarks {u, d - quarks) with current masses m„ ~ ~ 7 MeV acquire 
the dynamical masses fi^ — f^d — 320 MeV. This phenomenon is important 
for light pseudoscalar mesons as they possess the Nambu-Goldstone nature. 
To get low masses of pseudoscalar mesons one needs to take into account spin 
interactions of quarks. In [8,9] CSB was explained by the nonvanishing den- 
sity of quark (quasi) zero modes in the framework of QCD. Then the familiar 
PCAC (partial conservation of axial vector current) theorems and the soft 
pion technique are reproduced. There are different stochastic vacuum con- 
figuration which are responsible for CSB: instantons (anti-instantons), pieces 
of (anti-) self-dual fields (for example, torons or randomly distributed lumps 
of field), and others. The necessary requirement is to have zero fermion 
modes. The condensation of zero modes leads to CSB. The confinement of 
quarks does not allow them to be observed; i.e. quarks cannot move outside 
of hadrons in large distances relative to each other. This was confirmed by 
Monte-Carlo simulations and experiments. Both nonperturbative effects of 
strong interactions can be explained by introducing stochastic gluon vacuum 
fields with definite fundamental correlators [6,7]. Then the linear potential 
between quarks appears and it provides the confinement of quarks. Besides, 
Regge trajectories are asymptotically linear with a universal slope [7]. So the 
method of vacuum correlators in nonperturbative QCD and the dynamics of 
zero modes give the explanation of the double nature of light pseudoscalar 
particles (pions, kaons and others) as Nambu-Goldstonc particles and as the 
quark- ant iquark system with the confining linear potential. It should be no- 
ticed that confinement prevents the delocalization of zero modes over the 
whole volume [8], i.e., stabilizes the phenomena of CSB. 

In the present approach we make some assumptions. So we treats spin 
degrees of freedom as a perturbation and therefore it is questionable to apply 
this scheme to pions and kaons. Only p and K* mesons are considered here 
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because the energy shift for them due to the hyperfine spin interaction is 
below 100 MeV [4]. Here short-range spin-orbital L ■ S and spin-spin Si • S2 
interactions are not taken into account. Besides we neglect the Coulomb- 
like short-range contribution due to the asymptotic freedom of QCD. This 
contribution is important only for heavy quarkonia [5]. We imply also that 
in the presence of light quarks the structure of the vacuum yields an area law 
of the Wilson loop integral. The restriction of the leading Regge trajectories 
is used as we consider here only p and K* mesons. 

It is important to calculate different intrinsic characteristics of hadrons 
on the basis of QCD string theory and to compare them with experimental 
values. It will be the test of this scheme. The charge radius (and elec- 
tromagnetic form-factors) and electromagnetic polarizabilities of mesons are 
fundamental constants which characterize the complex structure of particles. 
These values for some mesons are known from the experimental data and 
therefore the estimation of them is reasonable. So we can check our notion 
about the vacuum structure by coinciding experimental data and theoretical 
predictions. In [12] we made the crude estimate of the charge radius and the 
electric polarizability of mesons and nucleons. Here we evaluate the mean- 
squared radius (using the virial theorem) and the magnetic polarizability of 
p and K* mesons. 

The electromagnetic polarizabilities of hadrons a, (5 enter the induced 
electric D = aE and magnetic M = /9H dipole moments, where E, H are 
the strengths of electromagnetic fields. As a result there is a contribution to 
the polarization potential [13,14] as follows 

C/(«,/3) = -iaE^-l/3H^ (1.1) 

Electromagnetic polarizabilities are fundamental low-energy characteris- 
tics of strong hadron interactions and therefore they can be calculated in 
the framework of non-perturbative quantum chromo dynamics — QCD string 
theory. 

This paper is organized as follows. In Sec. 2 after describing the general 
background we derive the effective action for mesons in external electromag- 
netic fields. The ground state and charge radii of particles are found on the 
basis of exact solutions and the virial theorem in Sec. 3. Section 4 contains 
the evaluation of the diamagnetic polarizabilities of p and K* mesons using 
the perturbative expansion in a small magnetic field. In the conclusion we 
made a comparison of our results with other approaches. 
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Units are chosen such that h — c — 1. 



2 Effective Action for Light Mesons 



To get the effective action for mesons in an external electromagnetic field we 
use the Green function of the quarks possessing spins in Minkowski space 
[see Eq. (A8) in the Appendix)]: 



S{x,y) ^ij^ ds J Dz [mi - -7^^^ {t)j Pexpji [^^2(0 - 



+eiZ,{t)A'^{z) + E^, (eiF;^ + (?F^,)]cit|$(a;, y), (2.1) 

where -2^(t) is the path of the quark with the boundary conditions ^^^(0) = 
y^, Zf,{s) = and 



$ (x, y) = P exp^ig j A^dz^ 



(2.2) 



is the path-ordered product (the parallel transporter), is the gluonic 
field and g is the coupling constant. Neglecting quark-antiquark vacuum 
loops and omitting the annihilation graph, the Green function of mesons 
(the quark-antiquark system) takes the form [7] 



G{x, x\ y, y) = tr {-t-,S{x, y)^ {y, y) j^Siy, x)^ {x, x)) , 



(2.3) 



were the brackets (...) are the averaging over the external vacuum gluonic 
fields with the standard measure exp [iS" (A)]. Inserting (2.1) into (2.3) we 
find the expression 

G{x, x; y, y) = tr / ds ds Dz(mi + --f^^Z/^ {t) 

Jo Jo Jz{0)=y V 2 



X [^^J^ Dz (^m2 - ^Tm^m (*)) expj^ [^i^^i) - + eii^(t)A^'(z) 



+eiE^,F„t iz^ 



tj-t^-^ 1x1/ 



dt + ij^ [l^lii) -ml + e2l^(f)<(^) + e^E^^F^l (z) 



dt 



X ( exp < igT, 



£ F^, {z) dt - F^, {z) dt^ I W{C)^ , (2.4) 



where Ps is the ordering operator of the spin matrices E^,^; ci, 62 are 
charges and mi, 7712 are current masses of the quark and antiquark; z^{t), 
z^{t) are the paths of the quark and antiquark with the boundary conditions 
z^{Q) = Vi^, Zi,{s) = x^, Zi^{0) = y^, Zf,{s) = and z^{t) = dz^/dt. Here we 
used the properties of 7 matrices: {75,7//} = 0, [2^,^,75] = 0. As compared 
with [7,17] we added the interaction of charged quarks with the external 
electromagnetic fields. The gage - and Lorenz-invariant Wilson loop operator 
is given by 



where Nc is the color number, and C is the closed contour of lines xx and 
yy connected by paths z{t), z{t) of the quark and antiquark. The Wilson 
operator (2.5) contains both the perturbative and non-perturbative interac- 
tions between quarks via gluonic fields A^. In accordance with the approach 
[7], spin interactions can be treated as perturbations. It is justified for p and 
K* mesons. To construct the expressions in spin interactions we write the 
relationship [7] 



where Sa/^i, {t) is the surface around the point 2;^ (t). The zeroth order 
in spin-orbit and spin-spin interactions corresponds to neglecting the terms 
gT,n^Fni, in Eq. (2.4). We suppose that mesons consist of quarks which move 
slowly with respect to the time fiuctuations of the gluonic fields (Tg ^ Tg). 
It is the potential regime of the string. Voloshin [19] and Leutwyler [20] 
remarked that in another case (Tg <^ Tg) the dynamics is nonpotential and 
the QCD sum rules can be used. We consider the case when the distance 
between quarks r » Tg. Monte-Carlo calculations [21,22] gave Tg ~ 0.2 -=-0.3 
fm. So we imply that the characteristic quark relative distance is r ^ 1 fm. 
This assumption will be confirmed below by the calculation of the quark- 
antiquark relative coordinate. 

The average Wilson integral (2.5) at large distances in accordance with 
the area law can be represented in Minkowski space as 




(2.5) 




(2.6) 



(W^(C)) =exp(zao5). 



(2.7) 
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where (Tq is the string tension and S is the area of the minimal surface 
inside of the contour C. The surface S can be parametrized by the Nambu- 
Goto form [23,24] 

S = £dT d(3^{w,w'^y-wlw% (2.8) 

where li;^ = duUfj^/dr, w'^ — dw/^/dp. Using the approximation [7] that 
the coordinates of the string world surface w^(t, P) can be taken as straight 
lines for the minimal surface we write 

w,{t, (5) = z,{t)(3 + z,{t){1 - P), (2.9) 

where r is implied to be the proper time parameter for both trajectories 
T — {tT)/s = {tT)/s. For uniform static external electromagnetic fields we 
have the representation of the vector potential through the strength tensor 

pel 

A:\z) = ^F^^lz,, Al^iz) = ^F^lz,. (2.10) 

The paths z^, are expressed via the center of mass coordinate and 
the relative coordinate [7], 

z^{T)^R^--^r^, z^{T)^R^ + -^r^ (2.11) 

o I o o I o 

with the boundary conditions for R^{t), r^{T): 

p /n^ /^l^M + /^2yM D (rr\ /^l^M + /^2^M 

RixK^) = : , RiJi^J- ) = : , 

The integration with respect to 2;^, in (2.4) is replaced by an integration 
over new variables i?^, r^. As r is a common time for the quark and antiquark 
(the time of the meson) the parametrization = (r, z), = (r, z) is possible 
[7]. This leads to the constraints: Ro{t) = r , ro(T) = 0. In accordance with 
the approach in [7] we introduce the dynamical masses //i, ^2 by relationships 
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Replacing the integration with respect to s, s in Eq. (2.4) by the integra- 
tion over d^i and dfj,2 with the help of (2.7) - (2.11) we find [12] the two-point 
function in zeroth order in the spin interactions 



G{x,x;y,y) = -T 
with the effective action 



d^i f°° dfi2 



2//i Jo 2/^2 



DRDr exp {iSe/f} 



^eff 



f 

Jo 



dr 



^2 1 / ^ 62 , 1~.^2 



2/ii 2p2 2 



J 



where 



(2.13) 



(2.14) 



Wf,^ Rij,+ [(5 - jJii/ini + A*2)]r^; A* = A*iA*2/(Aii + A*2) 

is the reduced mass of the quark-antiquark system, e = ei -|- 62, q — ei — e^.- 
As a first step we are interested here in the spinless part and therefore the 
preexponential terms {rrii + |7^i;^ (t)^, {rn2 — ITm^m (^) ^'^'^ constant 
matrix was omitted in (2.13). The expression (2.14) defines the 

effective Lagrangian for light mesons in external uniform static electromag- 
netic fields in accordance with the formula Seff = Jq Ceffdr. The expression 
(2.14) looks like a nonrelativistic one at F^l = 0, but it is not. The author 
of [7] showed that the relativism is contained here due to the jl dependence 
and the spectrum is similar to that of the relativistic quark model. 

The mass of the lowest states can be found on the basis of the relationship 

[25] 

J DRDrexp{iSeff} 

= (fi = ^^i^. r = X - S I exp{-iTM^, I = ^^14M, , = , _ p\ . 

\ A*i + Ai2 A*i + A*2 / 

(2.15) 

where the mass Al(/Xi, is the eigenfunction of the Hamiltonian. After 
that the Green function (2.13) is derived by integrating Eq. (2.14) over 
the dynamical masses /ii, ^2- In accordance with [7] we estimate the last 
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integration on d^i, djj,2 using the steepest descent method which gives a good 
accuracy when the Minkowski time T ^ oo. To have the correct formulas, 
it is necessary to go into Euchdean space and return into Minkowski space 
on completing the functional integration. We use this procedure. 



3 Ground State and Charge Radii 

The last term in (2.14) can be approximated with the accuracy of ~ 5 % [7] 
by the relation 

^1 



L 



^ dp. 



r2-(/5 (rxr)2~AA^. (3.1) 



It is the potential regime at low orbital excitations of the string when the 
orbital quantum number / is small. As the equalities Rq{t) = r, ro(T) = are 
valid, only three-dimensional quantities are dynamical. Taking into account 
(3.1), from Eq. (2.14) using the standard procedure we find the canonical 
three-momenta corresponding to the center of mass coordinate i?^ and the 
relative coordinate 

n. = ^ = (/.I + ^i2)Rk + If:Ir^ + iF^ir., 

TT, = ^ = Hrl + '-F:ir^ + ^-F^lR^. (3.2) 
The Hamiltonian 

Ti = -KkT-k + TlkRk — Ceff 

found from Eq. (2.14) with the help of Eqs. (3.1), (3.2) takes the form 

7i = ^ + ^ + ^^i±i^ + ^^i±^i?^ + ^/f. - ^(ER) 
2/xi 2/^2 2 2 2 ^ 2^ ^ 

-|(Er) + aoV^ (3.3) 
so that the equation for the eigenvalues is given by 

7i* = M(/xi,/X2)$. (3.4) 
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The terms containing the strength of the electric field in Eq. (3.3) describe 
the interaction of the dipole electric moment d with an external electric field. 
Using the definitions we have 

|(ER) + |(Er) = i(eiri + e2r2)E = dE (3.5) 

and the interaction energy of the electric dipole moment with a uniform 
static electric field is U = — dE. Bellow we investigate the case of a pure 
magnetic field when E = 0. The case E 7^ 0, H = was considered in 
[12]. Using Eq. (3.2) the equation for the eigenfunction $ of the auxiliary 
"Hamiltonian" 

n = n- ml/2iii - ml/2ii2 - {l^i + 1^2) /2 

is given by 

— f7r--(rxH)-^(RxH)V 

+— ^ r fn - ^(R X H) - ^(r x H)V + aoV^]^ 

= e(/i,H)$, (3.6) 

where e(/x, H) is the eigenvalue. In accordance with the Noether theorem 
we come to the conclusion that the canonical momentum 11 corresponding 
to the center of mass coordinate is a constant, i.e., 11 = const. Therefore it 
is possible to choose the condition 11 = 0. Putting 11 = and R = into 
Eq. (3.6) we arrive to the equation 

= efcH)<l>. (3.7) 

The second term in (3.7) describes the effect of the recoil of the string. 
Such a term appears also in non-relativistic models [26,13,14]. If we put 
R = in Eq. (2.13), this term would not appear [12]. 

In quantum theory instead of the path integration in r we can use the 
replacement TTfe — > —id/drk- We can apply Eq. (3.7) to the leading trajec- 
tories with light quarks with masses mi = m2 = m, //i = //2 = A* (/i = A*/2) 
for p meson and when mi 7^ m2, i^i 7^ /i2 for K* meson. 
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An external magnetic field splits the energy levels like the Zeeman effect 
for atoms. The difference with our case is we describe here the light quark- 
antiquark system in center mass system (c.m.s.) with the linear potential 
between quarks. Therefore the spectrum of the energy has other levels. 

We can consider a small external magnetic field so that here perturbative 
theory can be applied. We receive a first approximation when the external 
field H is switched off (H = 0) and the equation for the eigenvalue is given 
by 

(4^ + "°^)* = ^<''^*- 

Equation (3.8) gives the discrete values of the energy e(/i) due to the 
shape of the potential energy. The numerical solution of equation (3.8) was 
obtained in [27]. It is useful to find the solution to Eq. (3.8) for the ground 
state in an analytical form. After introducing the variables pk — {2flaoy^^rk, 
e{jl) - {2jl)-^'^al''^a{n) [7], Eq. (3.8) becomes 

+ pU(p)=a(n)$(p). (3.9) 



V 9p'. 

The solution to Eq. (3.9) may be chosen in the form $(p) = R{p)Yim{0, 0), 
where Yi^^O, (p) are spherical functions. After setting the variable R{p) = x/ P 
we come to the equation for the radial function 

X"(P) + («H - P - ^^^) Xip) = 0, (3.10) 

where x" (p) = d'^xip)/9p'^, and / is an orbital quantum number. The 
solutions to Eq. (3.10) for the ground state / = are the Airy functions 
Ai{p — a(n)), Bi{p — a{n)) [28]. The finite solution to Eq. (3.10) at p — > oo 
{I = 0) is 

x{p)^NAi{p-a{n)). (3.11) 
The constant N can be found from the normalization condition 

X\p)dp=l. (3.12) 

The requirement that this solution satisfies the condition 

X(0) = NAi{-a{n)) = 
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gives the Airy function zeroes [28] , 



a(l) = ai = 2.3381, a(2) = 02 = 4.0879 

and so on. The principal quantum number n = rir + l + l, where is the 
radial quantum number which defines the number of zeroes of the function 
x(p) at p > 0. For the ground state we should take the solution (3.11) at 
a{n) = oi (here = 0, / = 0): 

Xo{p)^NoAi{p-ai). (3.13) 

Now let us estimate the mean-squared radius for the state which is de- 
scribed by the function $ [the solution of equation (3.8)]. Multiplying Eq. 
(3.8) by the conjugated function and integrating over the volume we find 
the relations 

(T) + {U) = 6(/i), 
(T) = ^*dl^dV, (U) ^aoj Vr^^*^dV. (3.14) 

It is seen from Eqs. (3.14) that the mean potential energy {U) = ctq (^y/r^ 

is connected to the mean diameter (^a/t^^ (because r is the relative coordi- 
nate and quarks move around their center mass), which defines the size of 
mesons. In accordance with the virial theorem [29] we have the connection 
of the mean kinetic energy with the mean potential energy 

2{T) = k{U), (3.15) 

where k is defined from the equality U{Xr) — A*^C/(r). In our case of the 
linear potential k — 1 and from Eqs. (3.14), (3.15) we get 

{U) = ^e{fl) = ^{2fl)-'/'a'J\{n). (3.16) 

The use of the steepest descent method for the estimation of the integra- 
tion in fjL (at H = 0) leads to the conditions [7] : 

dMifxi,H2) ^ Q dM{fj.i,H2) ^ Q .3 
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where the mass of the ground state M.{fJ,i, 112) is given by [see Eqs. (3.3), 
(3.4)] 

/^2) = g + g + + (2^)-/3,„^/3,(,). (3.18) 

Here we consider the more general case as compared with [7] when /^i 7^ 
A''2 {nil 7^ T(i2). This case is reahzed for K* mesons. It is assumed that the 
current mass of u,d quarks (m„ = 5.6 ± 1.1 MeV, = 9.9 ± 1.1 MeV [30]), 
mi is much less than the dynamical mass /ii {^i ~ 330 MeV), i.e., mi -C /^i 
and the mass of s quarks 1112 (m^ = 199 ± 33 MeV [30]), is comparable with 
Hi but 1712 < A*i- Using these assumptions we neglect the term mi/2/xi in 
Eq. (3.18) and from Eqs. (3.17) have the equations 

(2/i(7o)^/^a(n) = Smj + {2flaof/^a{n) = (3-19) 

From Eqs. (3.19) we arrive to the expression for the dynamical mass 112 
(for s quarks): 

pi2 = \/pii + ml (3.20) 

To find III the perturbation in the parameter m\/ ix\ will be assumed. 
Using the relation 112 — A*i(l + to^/ (2//i)) which is obtained from Eq. (3.20) 
and the definition of the reduced mass jl — //1//2/ (a*i + A*2) from Eqs. (3.19) 
we arrive at the equation 




(3.21) 



In the zeroth order we come to the value //q = A*i = V^('^(^)/3)^^^ \J]- 
The next order gives the relationship 

-V.(t^)^-|(^)'''. (3..) 

In a particular case m2 = we arrive at fii = fi2 = fJ^o = 
[7]. The value of the string tension (Tq = 0.15 GeV^ was found from a 
comparison of the experimental slope of the linear Regge trajectories a' — 
0.85 GeV~^, and the variable a' = l/8o"o [7]. It leads for the lowest state 

= 0, / = 0, a(l) = 2.3381 to the value /xq = 321 MeV [7]. This means that 
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for p mesons when mi = m2 = we have the dynamical masses oiu, d 
quarks Pi = fi2 = /^o- For K* mesons using Eq. (3.20) and m2 = ~ 200 
MeV [30] from Eq. (3.22) we get the reasonable values 



/ii - 337 MeV, /la ~ 392 MeV, fi ~ 181 MeV. 



(3.23) 



Inserting the equation {U) = yo"o v r^y into the left-hand side of Eq. (3.16) 
one gives the expression 



o(2/U(7-o 



1-1/3 



an. 



(3.24) 



From Eqs. (3.20), (3.22) using the first order in the parameter rn^/ p\ we 



find 



( 3m2^ 



an 



3/4' 



(3.25) 



Equation (3.24) with the help of Eq. (3.25) gives an approximate relation 
for the mean relative coordinate: 



' a{n) 



3/4 



3m2 / 3 



3/2' 



-1/3 



(3.26) 



For p mesons putting nia = in Eq. (3.26) we arrive at 



an 



3/4 



(3.27) 



The same expression (3.27) was found in [12] using another method. With 
the help of the definition of the center of mass coordinate we can write an 
approximate relation for the mean charge radius of p-mesons: [] 



PI 



(3.28) 



At (To = 0.15 GeV^ [7,11] and a(l) = 2.3381 Eqs.(3.27), (3.28) give 



0.42 fm 



(( 



0.84 fm 



(3.29) 



1 The relationship (r^) ~ (vr^) is confirmed by the numerical calculations. 
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The value (3.29) characterizes the radius of the sphere where the wave 
function of the p meson is concentrated (remember that r is the distance 
between quarks). We know only the experimental data for tt^ mesons which 
have the same quark structure as mesons: 

rl^)^^^ = (0.44 ± 0.02) fm^ (^/(^,,, ^ 0.66 fm) [31]. 

For calculating the relative coordinate of K* mesons we should use Eq. 
(3.24) or (3.26) with the conditions /ii = /i„ and p2 = Eqs. (3.23). As a 
result formula (3.24) gives the value of the mean relative coordinate of K* 
mesons: 

'Vr2\ ^ = 0.79 fm. (3.30) 



/ K* 

With the help of this relation we can estimate the mean charge radius of 



K* mesons: 



The experimental data of the mean charge radius of mesons having 
an analogous quark structure as K* mesons are 



r^±) = (0.53 ± 0.05) fm [32], 
r|±) = (0.34 ±0.05) fm^ [31] 



and for neutral mesons 

^J{rlo) = (0.28 ± 0.09)fm [32] . 

So expression (3.26) gives a reasonable value for the charge radius of K* 
mesons, although we need the experimental data of the charge radius of K* 
mesons . 

The first perturbativc onc-gluon exchange contribution to Hamiltonian 
(3.3) determines the spin-spin correction such as the Breit-Fermi hyperfine 
interaction [9]. The spin-spin interaction is important to explain the Nambu- 
Goldstone phenomenon which takes place for the I = s = channel. 
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4 Perturbative Expansion and Magnetic Po- 
larizabilities 

To calculate the magnetic polarizability of mesons in accordance with Eq. 
(1.1) we should know the Hamiltonian depending on the magnetic field H. 
Prom Eq. (3.7) we arrive at the expression of the auxiliary "Hamiltonian": 

n = - — ^ + — HL +(^ + fr^H^ - (rH)^ 

+aoVr^, (4.1) 

where L = —i{rxd) {8^ = d/drk) is the angular momentum. Considering 
the external magnetic field H = (0, 0, H) expression (4.1) is rewritten as 

^ = -^1^2 + + hp + ^ (rl + rt) + <7oV^, 4.2 

2/i drf SjjL \AiJi Hi + IJL2 J 2,2 ^ ' 

where L3 = i{rid\ — rxd'i) is the third projection of the angular momen- 
tum. 

Now let us consider the magnetic polarizability of mesons on the basis of 
perturbative theory. We can rewrite Eq. (4.2) in the form 

n^no + —Ls + -r + :^ , 4.3 

Sfj, \4:jj, 111 + 1^2 J 32 

where 9 is the angle between H and the relative coordinate r and the free 
Hamiltonian is given by 

For the small magnetic fields H the second and third terms of Eq. (4.3) 
can be considered as a perturbation. Then using the standard perturbative 
method [33] , we find a shift of the energy in the state | n) with the accuracy 
of second order in H: 

.r. / , ei^, fe^ q^ \ i/Vsin^^ , , 

Sji \4/x III + 1^2 J 32 
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, — £/ / 

n " 

For the ground state I — {s state) the first and third terms of Eq. (4.4) 
do not give a contribution to the energy because L3 | 0) = 0. Taking the 
mean value and using the condition (l/47r) / sin^ OdQ = 2/3, from Eq. (4.4) 
we come to 

V) . (4.5) 




II1 + II2J 48 

Comparing Eq. (4.5) with Eq. (1.1) we find the magnetic polarizabihty 
of hght mesons 

0=~(^^^)(^')- (4.6) 



24 \4/i jJ^i + jJ,2, 

It should be noticed that here (r^) means the mean-squared relative co- 
ordinate of the quark- antiquark system. Expression (4.6) is similar to the 
Langevin formula for the magnetic susceptibility of atoms. It is seen that 
we have here only the diamagnetic polarizability as /? < 0. To calculate the 
paramagnetic polarizability one needs to take into account the interaction 
of the meson spin with the magnetic field. Using the approximate relation 

(r^) ~ (V^"^ , parameters e = Ci + 62, = Ci — 62, expression (3.27), and 
the dynamical masses of m, d quarks /xi = /X2, fl = n/2we find from Eq. (4.6) 
the relationship for the diamagnetic polarizability of p mesons: 

Calculating Eq. (4.7) for charged p mesons at (Jq = 0.15 GcV^, p = ji^ = 
321 GeV, ei = 2e/3, 62 = e/3, = 1/137, n = 1 in Gaussian units one takes 
the value 

(5p ~ -0.8 X 10-^ fm^ (4.8) 
In rationalized units, the polarizability is 47r times greater. 
For K* mesons, Eq. (4.6) at the values (3.23), (3.31) (^(r^) ~ (^)^) 
leads to the magnitude 

(3k* = -0.57 X 10"^ fm^ (4.9) 

Unfortunately there are not experimental data of p and K* meson polar- 
izabilities yet. 
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5 Conclusion 



The QCD string theory allows us to estimate the mean squared radii and 
magnetic polarizabilities of p, K* mesons. These quantities were derived as 
functions of the string tension which is a fundamental variable in this ap- 
proach. It is not difficult to calculate the magnetic polarizabilities of excited 
states of mesons using this approach. For that we should take the quantum 
numbers = 1, / = {n = 2) and evaluate the mean relative coordi- 
nate in accordance with Eq. (3.26). Then Eq. (4.6) gives the necessary 
polarizabihties. To have more precise values of the meson electromagnetic 
characteristics one needs to take into account spin corrections. Especially 
it is important for light pseudoscalar mesons (vr, K mesons). In principle it 
is possible to receive spin-orbit and spin-spin interactions using the general 
expressions (2.4), (2.6). 

The Nambu-Jona-Lasinio (NJL) model [34] having a good basis in the 
framework of QCD [35] describes chiral symmetry breaking but not the con- 
finement of quarks [7,36]. Besides this model has free parameters and the 
calculated polarizabilities of mesons [37] are parameter dependent. 

The instanton vacuum theory (IVT) developed in [39-40] does not give the 
confinement of quarks phenomenon [7]. This theory is like the NJL model 
[36] and takes into account only chiral symmetry breaking. Therefore the 
calculation of the meson electromagnetic polarizabilities on the basis of the 
IVT gave the similar result [41] as in the NJL model. 

All this shows that the theoretical evaluation of the charge radii and the 
magnetic polarizabihties of p, K* mesons is possible on the basis of a good 
description of chiral symmetry breaking and the confinement of quarks in the 
framework of QCD string theory but with some approximations and model 
assumptions. Naturally that theory was derived using the nonperturbative 
QCD, i.e., first principals of QCD. 

APPENDIX 

In this appendix we derive the one-quark Green function using the Fock- 
Schwinger method. Starting with the approach [15] and introducing external 
electromagnetic and ghionic fields we write the Green function of quarks 
which possess spins in Minkowski space 

s{x, y) = (x \ {d + ^i) \ y) = \ (^1 ~ -d) ~ -D^) I y) 
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(Al) 



where ei and mi are the charge and mass of the quark, D = 'y^D^, 
= dn — ieiA^J^ — igA^, A^ = A^A"; and A" are the Dirac and Gell- 
Mann matrices, respectively; A^J; and A^ are the electromagnetic and gluonic 

vector potentials, respectively. The inverse operator (ml — D'^^ can be 
represented in the proper time s [16]: 



{A2) 



Using the properties of Dirac matrices {7,^, 7^} = 26^, we find the squared 
operator 

D'^Dl + E^, (eiF;^ + gF,,) , (A3) 

where 

F^y = dfj,A^ - duAf, - ig[A^, A^], 



are the spin matrices, and F^j^, F^i, arc the strength of electromagnetic 
and gluonic fields, respectively. Inserting relationship (A2) into Eq. (Al) 
with the help of Eq. (A3) we get 



S{x^ y) — i J ds{x \ (rrii — exp 



mi 



Di 



(A4) 



The exponent in Eq. (A4) plays the role of the evolution operator which 
defines the dynamics of the "Hamiltonian" m\ — D'j^ — E^,^ {^iF^l + gF^iy) 
with initial | y) and final {x \ states where s means the proper time. Therefore 
it is convenient to represent the matrix element in Eq. (A4) as a path integral 
[16]: 



S{x,y) 



I / dsN / DpDzP (mi - d) exp i 

Jo Jz(0)=y ^ ^ L Jo 



fz(s)=x 



dt 



(p^ - eiA"^ - gA^y + E^, (eiF^^ + gF^,) 



{A5) 
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where (t) is the path of a quark with the boundary conditions z{0) — y, 
z{s) = X, D = (pfj^ — eiA^ — gA^^ and P means ordering; A?" is a constant 
which is connected to the measure definition and it will be chosen later. The 
path integration over the momenta can be rewritten in the form (see [7,18] 

N J Dp (mi - exp |i ^ dt p^z^ - {p^ - eiA"^ - gA^f | 



J Dp exp i dt{p^z^) 
X exp |-i ^ dt [p^ - eiA'^ - gA^f^ 
N J Dp (mi - ^-ff,z^^ exp |z ^ dt p^z^ - (p^ - e^A^^ - gA^^j^ | 



j Dp 



mi - -^i^Zf, J exp <; z 



dt 



-pl + lil + {(^iAi + gA^)z, 



(AQ) 

In Eq. (A6) we used the integration by parts (see [18]) and made a 
continuity of shifts p^i ^ p,j. + ^lA^j + gA^ and then p^ ^ p^ + 2^/2. The 
constant N in Eq. (A6) is defined hy the relation 



N I i^pexpj-i = 



{A7) 



Taking into account Eqs. (A6), (A7) we find from Eq. (A5) the Green 
function of the quark 



S{x,y) 



ds 



z(s)=x 
ziO)=y 



Dz [mi - --ft^Zn (t) P exp 



dt 



+ (eiA'l + gA^ z^ (t) + E^, (eiF^l + gF^,) |. 



-zl (t) - m\ 
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